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ABSTRACT 

The radii of some transiting extrasolar giant planets are larger than would be expected by the stan- 
dard theory. We address this puzzle with the model of coupled radius-orbit tidal evolution developed 
by Ibgui & Burrows (2009). The planetary radius is evolved self-consistently with orbital parameters, 
under the influence of tidal torques and tidal dissipation in the interior of the planet. A general feature 
of this model, which we have previously demonstrated in the generic case, is that a possible transient 
inflation of the planetary radius can temporarily interrupt its standard monotonic shrinking and can 
lead to the inflated radii that we observe. In particular, a bloated planet with even a circular orbit 
may still be inflated due to an earlier episode of tidal heating. We have modified our model to include 
an orbital period dependence of the tidal dissipation factor in the star, Q'„ oc P 7 , — 1 ^ 7 ^ 1. With 
this model, we search, for a tidally heated planet, orbital and radius evolutionary tracks that fall 
within the observational limits of the radius, the semimajor axis, and the eccentricity of the planet in 
its current estimated age range. We find that, for some inflated planets (WASP-6b and WASP-15b), 
there are such tracks; for another (TrES-4), there are none; and for still others (WASP-4b and WASP- 
12b), there are such tracks, but our model might imply that we are observing the planets at a special 
time. Finally, we stress that there is a two to three order-of-magnitude timescale uncertainty of the 
inspiraling phase of the planet into its host star, arising from uncertainties in the tidal dissipation 
factor in the star 

Subject headings: planetary systems — planets and satellites: general — planets and satellites: indi- 
vidual (WASP-4b, WASP-6b, WASP-12b, WASP-15b, TrES-4) 



1. INTRODUCTION 

The more than 60 transiting extrasolar giant planets 
(EGPs) discovered so far 1 offer a unique opportunity to 
test and improve the models of the structure and evo- 
lution of these bodies. The mass and radius of such 
planets can be inferred from a combination of radial ve- 
locity and transit lightcurve measurements that break 
the planet mass-inclination angle degeneracy. A large 
theoretical effort has been undertaken for more than a 
decade now to model and understand the evolution and 
the radii of transiting planets (Guillot et al. 1996; Bur- 
rows et al. 2000; Bodenheimer et al. 2001; Burrows et al. 
2003; Baraffe et al. 2003; Bodenheimer et al. 2003; Gu 
et al. 2003; Burrows et al. 2004; Fortney & Hubbard 2004; 
Baraffe et al. 2004; Chabrier et al. 2004; Laughlin et al. 
2005; Baraffe et al. 2005, 2006; Burrows et al. 2007; Fort- 
ney et al. 2007; Marley et al. 2007; Chabrier & Baraffe 
2007; Liu et al. 2008; Baraffe et al. 2008; Ibgui & Bur- 
rows 2009; Miller et al. 2009; Leconte et al. 2009; Ibgui 
et al. 2009). 

The radius of a gas giant planet depends on many phys- 
ical effects that are particular to a given planet-star sys- 
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tern, including the mass and age of the planet; the stellar 
irradiation flux and spectrum; the composition - in par- 
ticular, the heavy-element content - of the atmosphere, 
the envelope, and the core; the atmospheric circulation 
that couples the day and the night sides; and any pro- 
cesses that could generate an extra power source in the 
interior of the planet, such as tidal heating. Moreover, 
the transit radius effect (Burrows et al. 2003; Baraffe 
et al. 2003) has to be considered in order to infer the 
transit radius from the planet's physical radius. There- 
fore, a custom evolutionary calculation is the most ap- 
propriate way to determine a theoretical transit radius, 
in order to compare it with the observed transit radius. 

The objective of the present paper is to test the coupled 
radius-orbit tidal evolution model developed by Ibgui 
& Burrows (2009) on some recently discovered inflated 
planets. We present evolutionary tracks for WASP-4b 
(Wilson et al. 2008; Southworth et al. 2009; Gillon et al. 
2009b; Winn et al. 2009a) and WASP-12b (Hebb et al. 
2009). We have also tested the model on TrES-4 (Mandu- 
shev et al. 2007; Sozzetti et al. 2009), WASP-6b (Gillon 
et al. 2009a), and WASP- 15b (West et al. 2009). Some 
other inflated planets have been discovered recently, such 
as HAT-P-13b (Bakos et al. 2009), WASP-17b (Ander- 
son et al. 2009), and CoRoT-5b (Rauer et al. 2009). We 
might apply our formalism to them in the future. 

The idea of exploring tidal heating as an explanation 
for the inflated radii was originally formulated by Boden- 
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heimer et al. (2001). They suggested an excitation mech- 
anism to sustain a nonzero eccentricity, for example a 
planetary companion (Bodcnheimcr et al. 2003; Mardling 
2007). To date, two transiting EGPs are known to be 
accompanied by a companion, HAT-P-13b (Bakos et al. 
2009) and HAT-P-7b (Pal et al. 2008; Winn et al. 2009b). 
This reinforces the plausibility of such a scenario. 

Batygin et al. (2009) have coupled a three-body tidal 
orbital evolution model with a model of the interior struc- 
ture of HAT-P-13b. They found a quasi-stationary solu- 
tion and possibly consistent core masses, radii, and tidal 
heating rates. Assuming, as did Liu et al. (2008), that the 
systems are in a quasi-stationary state, Ibgui et al. (2009) 
provide, for each of the systems that they have stud- 
ied (WASP-6b, WASP-12b, WASP-15b, TrES-4, HAT-P- 
12b) and for each of the associated atmospheric opacities 
(solar, lOxsolar) a relation between the heavy-element 
core mass M core and the ratio e 2 /Q' p , where e is the or- 
bital eccentricity and Q' p the tidal dissipation factor in 
the planet. This constraint results from a degeneracy 
between the dissipation heating rate in the interior of 
the planet (which increases the radius) and the mass of 
a possible heavy-element central core (which shrinks the 
radius). 

For close-in EGPs (orbital separation < 0.1-0.15 AU), 
tidal torques are strong enough such that they can result 
in planetary orbital evolution, and they produce tidal 
heating (dissipation) inside the planet. Such tidal effects 
were first suggested for transiting EGPs by Jackson et al. 
(2008b, c,d). Jackson et al. included the tides raised on 
the star and the tides raised on the planet. They found 
tidal rates close to the levels that Burrows et al. (2007) 
proposed to maintain the observed radii of some transit- 
ing EGPs. Ibgui & Burrows (2009) described a model 
that couples the two consequences of these tidal effects 
- planetary radius evolution and orbit evolution. They 
tested their model on HD 209458b and found an expla- 
nation for the radius of this planet. Note that they also 
showed that a supersolar metallicity of the planetary at- 
mosphere, without invoking tides, can explain the ra- 
dius. Miller et al. (2009) applied a similar method to all 
the transiting EGPs, albeit with simplified models for 
the atmospheres of the planets and parent stars, and a 
restricted range of possibilities for the tidal dissipation 
factors Q' . We have chosen a more detailed approach 
by adopting customized atmospheric models and an ex- 
tended range for Q' . Therefore, due to the complexity of 
the atmospheric calculations, we have selected a couple 
of planets. We believe that both approaches are comple- 
mentary The one followed by Miller et al. (2009) pro- 
vides a global estimate of the possibilities for matching 
the observed radii of the transiting EGPs, while our ap- 
proach, more precise and therefore applied to fewer plan- 
etary systems, focuses on more specific issues such as the 
influence of the atmospheric opacity of the planet, a more 
detailed model for the tidal dissipation in the star, and a 
phenomcnological study of all the qualitatively possible 
behaviors of the evolutionary curves (Ibgui & Burrows 
2009). The application of our model to a subset of in- 
flated transiting EGPs is the subject of this paper. 

The paper outlines, in Section 2, the main assumptions 
of our coupled radius-orbit tidal evolution model, with a 



summary of the basic phenomenological results obtained 
from the previous study by Ibgui & Burrows (2009). 
It also explains our upgraded modeling of the tidal 
dissipation factor in the star. Section 3 demonstrates 
some additional generic results, such as the effect of 
M COTe , Q'p, and the initial semimajor axis. For each of 
the following planets, TrES-4, WASP-4b, WASP-6b, 
WASP-12b, and WASP-15b, we search for evolutionary 
tracks that fall within the observational limits of the 
radius, the semimajor axis, and the eccentricity of the 
planet in its current estimated age range. Our results 
are described in Section 4. We have not found any such 
track for TrES-4. We have found solutions for WASP-6b 
and WASP-15b. The cases of the planets WASP-4b and 
WASP-12b, for which we have coupled models that fit, 
are more interesting. Therefore, we present evolutionary 
curves for these planets. In fact, the solutions that we 
obtain for these two planets are valid only for very short 
age ranges in comparison with the estimated ages of the 
planets. This would imply that we are observing both 
planets at a very special time in their evolution, which 
would be a priori unlikely. In Section 5, we discuss the 
plunging timescale of a planet into its host star. Its 
uncertainty can span two to three orders of magnitude. 
We summarize our results in Section 6. 



2. THE COUPLED RADIUS-ORBIT EVOLUTIONARY 
MODEL AND THE PERIOD DEPENDENCE OF THE 
TIDAL DISSIPATION FACTOR Q' 

2.1. The Phenomenology of the Coupled Radius-Orbit 
Evolution 

The model that we employ assumes a two-body grav- 
itational and tidal interaction that consistently couples 
the evolution of the radius with the orbit of the planet. 
It includes the tides raised on the planet and the tides 
raised on the star, along with stellar irradiation and de- 
tailed model atmospheres. The planetary radius evolves 
as a result of the competing influences of tidal heating 
in its convective interior (due to the dissipation of or- 
bital energy) and radiative cooling from its surface. The 
most interesting phenomenological result that we have 
obtained is that, for strong enough tides, the planet's ra- 
dius can undergo a transient phase of inflation that tem- 
porarily interrupts its monotonic shrinking and resets its 
evolutionary clock (Ibgui & Burrows 2009). Moreover, 
we have demonstrated that, due to thermal inertia, an 
earlier episode of tidal heating can result in an inflated 
radius at the current age of the planet, even though its 
current orbit has nearly circularized. 

2.2. Formalism and Assumptions: a Summary 

The formalism, assumptions, and computational tech- 
niques have been extensively described in Ibgui & Bur- 
rows (2009). We summarize them here. 

The planetary structure consists of a gaseous (H 2 , He) 
iscntropic envelope (helium mass fraction Y = 0.25) de- 
scribed by the equation of state of Saumon et al. (1995). 
It may also contain an inner heavy-element core. The 
basic effect of the core is to shrink the planetary radius 
(Burrows et al. 2007; Liu et al. 2008; Ibgui et al. 2009). 
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However, in the context of coupled radius-orbit evolu- 
tion, its effect is more subtle, as will be discussed in §3. 
We restrict ourselves to a solar atmospheric opacity of 
the planet. The effects due to a higher opacity, namely 
an enhanced and accelerated transient phase of radius in- 
flation, have been described in Ibgui & Burrows (2009). 

The giant planet radius evolution is modeled with a 
Henyey code (Burrows et al. 1993, 1997). The radiative 
cooling from the surface of the planet is linked to bound- 
ary conditions (Burrows ct al. 2003). The latter incorpo- 
rate realistic irradiated planetary atmospheres calculated 
by CoolTLUSTY, a variant of TLUSTY (Hubeny 1988; 
Hubeny & Lanz 1995). The specific behaviors of such ir- 
radiated atmospheres are an active field of research. For 
example, an extra absorber in the upper-atmosphere may 
produce a temperature inversion (Hubeny et al. 2003; 
Burrows et al. 2008a, b; Fortney et al. 2008; Knutson 
et al. 2008; Spiegel et al. 2009c). We do not incorporate 
such an extra absorber in the models of this paper. High 
mctallicity EGPs, such as Neptune-mass planets (Spiegel 
et al. 2009a), have their own unique characteristics. We 
calculate a customized host star spectrum by interpola- 
tion of the Kurucz (1994) models at the actual effective 
temperature and gravity of the star. Such a customized 
approach can only improve the reliability of the model. 

The equations governing the tidal evolution of the or- 
bital eccentricity e, the semimajor axis a, and the tidal 
heating rate are the ones adopted by Ibgui & Burrows 
(2009). We rewrite them here for the sake of com- 
pleteness (Goldreich & Soter 1966; Kaula 1968; Peale 
& Cassen 1978; Murray & Dermott 1999; Bodcnheimer 
et al. 2001, 2003; Gu et al. 2004; Mardling 2007; Jackson 
et al. 2008c,d,a; Ferraz-Mello et al. 2008; Barnes et al. 
2009): 

1 de _ 
e dt 
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In the preceeding equations, G is the gravitational con- 
stant, M p , M*,R p , and i?» are the masses and radii of 
the planet and star, and Q' p and Q'^ are the tidal dis- 



sipation factors in the planet and in the star (Goldreich 
1963; Goldreich & Soter 1966). The planet radius is time- 
dependent, R p (t), but the star's radius is assumed to be 
constant. The principal assumptions underpinning these 
equations (more details can be found in Ibgui & Burrows 
2009) are as follows. We consider that this two-body in- 
teraction starts a few Myr after star formation, preclud- 
ing any interaction with a protoplanctary disk (Goldreich 
& Sari 2003) or with other planets (Ford et al. 2003; J uric 
& Tremaine 2008; Chatterjee et al. 2008; Ford & Rasio 
2008; Nagasawa et al. 2008). We also do not consider 
the Kozai interaction (Wu & Murray 2003; Wu 2003; Wu 
et al. 2007; Nagasawa et al. 2008). We assume that, after 
a few Myrs, the planet is close enough to its host star that 
tidal effects can be significant, typically using an "initial" 
semimajor axis up to 0.1 AU and we consider an "ini- 
tial" eccentricity a that can range from 0.0 to 0.8. We 
neglect stellar and planetary obliquities. We assume that 
the planet's spin is synchronized (tidally locked) with its 
orbital period, and that the star's spin rate is small com- 
pared with the orbital mean motion. These equations are 
developed to lowest-order in e (Goldreich & Soter 1966). 
They become less good approximations at higher values 
of eccentricity. We present results that include high val- 
ues of e for comparison with previous work (Jackson et al. 
2008c; Ibgui & Burrows 2009; Miller et al. 2009 use the 
same approximation). Higher-order terms may be con- 
sidered in the future (Mardling & Lin 2002). However, 
these equations rely on specific assumptions for the re- 
sponse of a body to tidal forcing. Tidal theory is still an 
active research field with many remaining ambiguities. 
Therefore, our approach is sufficient to describe the ba- 
sic phenomenology. These equations show that e and a 
can, given the above assumption, only decrease. 

Once the properties of both the star and the planet 
are fixed, the model is controlled by four free parameters 
(Q' P j Q'*,Si,ai) that we vary in order to fit simultaneously, 

' within the error bars of the age of the planet, the plane- 
tary radius, the eccentricity, and the semimajor axis. An 
important point to bear in mind is that the radius and 
orbit evolution are strongly nonlinear and depend sensi- 
tively on these four free parameters, e, is varied from 0.00 
to 0.80. en is varied from its current (measured) value to 
0.10 AU, the approximate distance beyond which tidal 

' effects are negligible in our cases. The tidal dissipation 
in the planet is ruled by Q' p , whose value is poorly con- 
strained. Experimental estimates provide 10 5 to 10 6 for 
Jupiter (Goldreich & Soter 1966; Yoder & Peale 1981), 
while theoretical arguments suggest around 10 5 to 10 7 
(Ogilvie & Lin 2004), and up to 10 8 for planets with cores 
(Goodman & Lackner 2009). Consequently, we have con- 
sidered the range 10 5 — 10 8 for Q' p . A particular model 
is specified for Q'„ as is explained in the next subsection. 



2.3. The Tidal Dissipation Factor in the Star 

The tides raised on the star play a specific role in that 
they are the major if not the only contributor to the 
planet's orbital evolution, as soon as the orbital eccen- 
tricity e is small enough or is zero. Moreover, for a zero 
eccentricity, the evolution of the semimajor axis is de- 
scribed analytically, as was first demonstrated by (Gol- 
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dreich 1963), as can also be directly derived from eq. (2): 

-i 2/13 



a 



117 G 1 ' 2 M p 

~^mF 2 q 



7(* 



to) 



(6) 



where ao is the semimajor axis at any time to after the 
eccentricity has become null. 

The subsequent evolution of the transiting EGPs from 
their observed current state is a puzzling issue. The sta- 
bility, but also the timescale, of their evolution are being 
investigated (Spiegel et al. 2009b; Levrard et al. 2009; 
Hellier et al. 2009; Hamilton 2009). It is well known 
that the tides raised on the star ultimately cause the 
planet to spiral into its host star and probably eventu- 
ally to be tidally disrupted (Rasio et al. 1996; Levrard 
et al. 2009; Jackson et al. 2009b,a; Ibgui & Burrows 2009; 
Miller et al. 2009). The associated timescale is higly de- 
pendent on Q' t . 

As an upgrade to the model used for our previous 
generic study (Ibgui & Burrows 2009), we consider that 
the tidal dissipation factor in the star, <3'„, can evolve 
with the orbital period of the system as follows: 



(7) 



where P is the orbital period of the system, Pq is a ref- 
erence orbital period (i.e. for which Q' t = 10 13 ). The 
exponent 7 is between —1 and +1. The theoretical mo- 
tivation for this range of 7 comes from the modeling of 
tidal dissipation in fluid bodies and the modeling of tur- 
bulent kinematic viscosity (Zahn 1966, 1989; Goldreich 
& Nicholson 1977). The range of 7 has been empiri- 
cally confirmed by Spiegel et al. (2009b) who show, based 
on a statistical study of the observed transiting planets' 
properties, that only such a range can lead to a station- 
ary rate of plunging planets into their host stars. The 
range for j3 is 5.0 - 8.0, resulting in Q'„ within 10 5 — 10 8 
at P — P . We adopt the same range as for Q' p . We 
choose Po — 10 days for the reference orbital period. In- 
deed, observational data from Meibom & Mathieu (2005) 
on solar-type binaries in the open cluster M35 suggest a 
value of <3'„ w 10 6 for a period of 10 days (Ogilvie & Lin 
2007). 

3. COUPLED EVOLUTION OF THE TIDALLY HEATED 
RADIUS AND THE ORBIT: THE EFFECTS OF M core , Q' p , 
AND a t 

Major generic features of the radius and orbital co- 
evolution of a close-in giant planet are described by Ibgui 
& Burrows (2009). Here, we present additional results: 
we evaluate the effect of M core , Q' p , m on the evolution 
of the radius and the orbit of a transiting EGP. Our 
statements can be demonstrated with equations (1,2,3). 
They have also been numerically tested on the "generic 
transiting system" employed by Ibgui & Burrows (2009), 
namely the HD 209458 system. Here, the tides raised on 
the star are neglected (<3* — > 00) for clarity's sake. Ev- 
erything else being equal, increasing the values of either 
M COTe , Q' p , or a,, results in delaying the appearance of 
the radius inflation peak and decreasing its value. It also 
results in faster circularization of the orbit, and the final 



semimajor axis a/ is reached faster. As for the latter, by 
virtue of conservation of the angular momentum of the 
system, a/ is the same whatever (for whatever M coro or 
Q'p), but a larger di results in a larger a/. 

The demonstration of the evolution of the radius in- 
flation peak and the orbit is based on the fact that a 
higher value of each of the three parameters (M coro , Q' p , 
di) results in a lower initial tidal heating rate (see eq. 3). 
It is straightforward for an enhanced Q' p or a*. Bur- 
rows et al. (2007) established that the larger the heavy- 
element core, the smaller the EGP radius. At the same 
time, the evolutions of e and a start with a lower rate. 
The tidal heating being initially lower, the ensuing tran- 
sient expansion phase is manifest to a lesser degree and 
later. 

The conclusion is that, when tidal heating is coupled 
with orbital evolution, the addition of a core does not 
necessarily result in a smaller radius. The result depends 
on the age of the system. In essence, at earlier age the 
planet with the larger core has the smaller radius, but 
it is the opposite at later age. 



4. COUPLED EVOLUTION OF THE TIDALLY HEATED 
RADIUS AND THE ORBIT: APPLICATIONS 

The validity of the model can be tested against the 
available data of the observed inflated transiting EGPs. 
Our objective is to find evolutionary tracks that fall 
within the observational limits of the radius, the semi- 
major axis, and the eccentricity of the planet in its cur- 
rent estimated age range. If we do find such tracks, we 
say that "we fit the planet" . The first application was 
presented in Ibgui & Burrows (2009) for HD 209458b. 
We were able to simultaneously fit the radius, the eccen- 
tricity, and the semimajor axis of this planet with the 
set (Q'p, Q'*,ei,ai)= (lO^ 55 , 10 7 , 0.77, 0.085 AU), where 
<3* is constant, and for a solar opacity. Note that if 
HD 209458b has 3 to 10 x solar opacity, we can explain 
its radius without invoking the tidal heating argument. 
In this paper, we present the results of our investigation 
for other inflated planets, WASP-4b (§4.2) and WASP- 
12b (§4.3). We assume a solar opacity, and no central 
heavy-element core. We adopt a Q'^ that varies accord- 
ing to eq. (7) with 7 = — 1 in order to smooth the ulti- 
mate plunging of the planet into its host star. We have 
also tested our coupled model for TrES-4, WASP-6b, and 
WASP-15b. Observational data are listed in Table 1 for 
the planets' properties and in Table 2 for the host stars' 
characteristics. 

For each planet, we have tested a large number of com- 
binations of the parameters (Q' p , Qi,ei, di). Recalling the 
Q'* evolutionary law given by eq. (7) in §2.3, the refer- 
ence orbital period is Po = 10 days. We took 7 values for 
the [3 parameter that determines Q' x at P = Po (f3 = 5.0 
to 8.0 in intervals of 0.5), 31 values for Q' p (log 10 (Qp) = 
5.0 to 8.0 in intervals of 0.1), 81 values for e (0.00 to 
0.80, in intervals of 0.01), and a certain number (39 for 
WASP-4b and WASP- 12b) of values of a (the observed 
value up to 0.10, in intervals of 0.002). This represents 
^680,000 evolutionary curves tested for each planet at a 
given opacity. Given the high nonlinearity of the evolu- 
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tionary equations, and the sensitive dependence on these 
parameters (see §2.2), this approach represents a fairly 
exhaustive exploration of all the possible combinations in 
order to select the ones that fit the observed parameters. 

In §4.1, we describe results of our models in the cases 
of WASP-6b, WASP-15b, and TrES-4. Perhaps the 
most interesting cases, however, are WASP-4b (§4.2 and 
Fig. 1) and WASP-12b (§4.3 and Fig. 2), which we can 
fit, assuming a solar opacity. 

4.1. WASP-6b, WASP-15b, TrES-4 

WASP-6b and WASP- 15b were discovered by the Wide 
Angle Search for Planets (WASP) survey (Gillon et al. 
2009a; West et al. 2009). The radius of WASP-6b is 
1.224±°;°^i?j, its mass is 0.503t^Mj, and its age is 
Gyr. The radius of WASP-15b is 1A28±° ° 7 7 7 7 Rj, 
its mass is 0.542i°;^°Mj, and its age is 3.9t?j Gyr. 
These two planets can easily be fit with little tidal heat- 
ing (log 10 (Qp) ^ 7.5) at solar opacity and for large age 
ranges. Moreover, the opacity effect is sometimes suffi- 
cient to explain some radii: WASP-6b can be fit with 
a 3 x solar opacity, without tidal heating (Ibgui et al. 
2009). Evolutionary tracks (not shown) that fit these 
planets' observed properties look very similar to those of 
HD 209458b, presented in Ibgui & Burrows (2009). 

TrES-4 was discovered by Sozzctti et al. (2009). Its 
radius is 1.783±3;g||J2j, its mass is 0.925t^Mj, and 

its age is 2.9±J;| Gyr. This planet is one of the most 
inflated transiting planets known, and its size is difficult 
to explain, even with tidal heating, at its estimated age. 
We were not able to simultaneously fit the radius, ec- 
centricity, and semimajor axis of TrES-4, at any of the 
opacities that we have tested (solar, 3 x solar, 10 x solar). 

4.2. WASP-4b 

WASP-4b was discovered by Wilson et al. (2008) and 
its parameters were further refined by Southworth et al. 
(2009), Gillon et al. (2009b), and Winn et al. (2009a). 
Its observed radius is 1.365+ 021 -Rj an d its mass is 
1.237tg;g^Mj for an age of 6.5±H Gyr. Its eccentric- 
ity is not well-constrained, with an upper limit of 0.096 
(Madhusudhan & Winn 2009). 

Examples of evolutionary curves that fit for WASP-4b 
are portrayed in Fig. 1. This figure consists of four pan- 
els that depict, versus the age (in Gyr) of the planet, 
the simultaneous evolution of its radius R p (Rj) (top left 
panel), its eccentricity e (top right), its semimajor axis 
a(AU) (bottom left), and log 10 ((2'„) (solid curves, left 
y axis) and the orbital period P(days) (dashed curves, 
right y axis) in the bottom right. Solar opacity is as- 
sumed for the planetary atmosphere. The radius evolu- 
tion without tides and a constant orbit is represented by 
a black curve in the top left panel. It shows that the stan- 
dard model prediction is quite far from the measurement. 
The difference is roughly 0.16 Rj, i.e. 12%. Among the 
tested <2'„ at P = P a (log 10 [Q'„(P )] = 5.0 to 8.0 in steps 
of 0.5), we find fitting solutions for log 10 [Q»(Po)] ^ 6.5. 
The smaller the the steeper the plunging slope of 



the semimajor axis. The results we show in Fig. 1 are 
among the ones that have the smoothest plunging orbits. 
Moreover, for all the fitting curves, log 10 (Qp) is greater 
than or equal to 6.9, which is a fairly high value. In the 
examples on Fig. 1, Q' p = 10 8 . Also, the initial eccen- 
tricity is large, ej = 0.80. Solutions with lower ej exist 
(the lowest is 0.45), but they are for Qi(P ) = 10 50 
or 10 5 5 , which arc the fastest plunging configurations. 
Three fitting evolutionary curves are plotted, for three 
different a,: 0.050, 0.052, 0.054 AU. We identify, when 
dj increases, the delay of the appearance, and the lower 
value, of the radius inflation peak, as stated in §3. The 
curves end with thick dots, where the periastron of the 
orbit reaches the Roche limit represented by a brown 
horizontal line in the bottom left panel. The ranges of 
ages for which simultaneous fits are obtained within the 
plotted la measurement limits for R p , e, a, are repre- 
sented by two orange vertical segments. These ranges 
are very narrow (< 0.2 Gyr) in comparison with the es- 
timated age of the planet, 6.5+ 2 ' 3 Gyr. It is awkward to 
suggest that we observe this planet at a very special tran- 
sitional period in its life. However, if we relax the fitting 
criterion (on R p , e, and a) from la to 2a or 3a, there 
would be a wider range of ages for which the evolution- 
ary tracks would fit the observed values. Therefore, our 
models would no longer imply that the present is such a 
special time in the life of this planet. The bottom right 
panel shows the evolution of Q'^, which increases with 
time, from log 10 ((3'„) ~ 6.8 to ~ 7.6 over the whole evo- 
lutionary path. This is because of its inverse dependence 
on the orbital period P, which decreases while the orbit 
circularizes, from ~4.5 days to ~0.8 days at the Roche 
limit. The fitting ranges are naturally at the measured 
period P m0 asurod ~ 1-3 days for which Q[ t k, 10 7 ' 4 . This 
increase of Q'„ corresponds to a decrease in the tidal dis- 
sipation in the star and to a smoother evolution of its or- 
bit. In the bottom left panel, we notice an apparent jump 
of the slope in the semimajor axis evolution, previously 
pointed out and justified by Miller et al. (2009) in their 
simulaton of the evolution of HD 209458b. This jump is 
respectively at 5, 6, 8 Gyr for a, t = 0.050, 0.052, 0.054 AU. 
It can be explained quasi-analytically using eq. (2) in a 
more compact form: 

a = a p e 2 + a* (l + ^e 2 ^j , (8) 

where a p and d* are the rates of evolution of a due to the 
tides, respectively raised on the planet and on the star 2 . 
Before this apparent jump, both tides contribute to the 
decrease of a, roughly by a comparable amount according 
to numerical tests. Then, during the rapid decrease of 
e as shown in top right panel of Fig. 1, the components 
of the sum (8) that depend on e fall even more rapidly 
because of the square dependence on e. Finally, after 
this short transitional phase, a evolves only as d», which 
is independent of e, leading to the analytical evolution 
described by eq. (6). 

Figure 1 shows that, once the orbit has circularized, 

2 Explicitly, the rates are defined by a p = —2K\ Vt n ;L, , and 

a / Q p 

■ !t - Rj 

a * ~ ~25 n2 P* 7 tt 72qT- 
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the radius still remains far above the value reached with 
no tides, as the top left panel demonstrates at the end 
of evolution depicted by the thick dots. This specificity, 
already mentionned by Ibgui & Burrows (2009), is due 
to thermal inertia. In sum, it is possible that an inflated 
planet with a circular orbit can be explained with tidal 
heating. 

4.3. WASP-12B 

WASP-12b, discovered by Hebb et al. (2009), is 
unique in many respects. It is the second-largest tran- 
siting planet to date, with a radius of 1.79^'qqRj. 
Its mass is IA1+ ']^Mj. Its estimated eccentricity is 

e = 0.049^0.015- ft i s tne m °st heavily irradiated 
transiting EGP with a flux at the substcllar point of 
9.098xl0 9 ergs cm~ 2 s _1 . It has one of the shortest 
orbital periods, P = 1.09142 days. Its orbit is very close 
to its Roche limit, 0.0221 AU, while its periastron, p = 
a(l - e), may be below this limit: 0.0207 < p < 0.0229 
(see Table 1). The planet is perhaps on the verge of be- 
ing tidally disrupted. It may also be losing mass due to 
Roche lobe overflow (Gu et al. 2003; Li et al. 2009). This 
phenomenon is not modeled here. 

In order to fit the radius, Miller et al. (2009) invoke 
a floor on the eccentricity and an extremely rapid ex- 
pansion of the radius when the planet starts to plunge. 
Our model provides fitting evolutionary curves without 
imposing a floor on e and without the rapid expansion. 
Figure 2, similar to Fig. 1, shows examples of fitting 
curves for a solar opacity atmosphere. As is done for 
WASP-4b, the radius evolution with no tides is drawn 
in black in the top left panel. The difference between 
the radii is ~0.50 Rj, that is 28%, much larger than for 
WASP-4b. Similar to WASP-4b, we obtain solutions for 
log 10 [Q'„(-Po)] ^ 6.5, and high values of Q' p are required, 
Q' p ^ lO 7 2 . We show in Fig. 2 three examples of fitting 
evolutionary curves with the smoothest plunging orbits 
(Q'*(Po) = lO 65 ), for three different a t : 0.053, 0.055, 
0.057 AU, and an initial eccentricity of ei = 0.73. Lower 
ei is possible, down to 0.53, but for lower Q'*{Pq) and, 
therefore, faster plunging orbits. The delay of the ap- 
pearance of the radius peak and its lower value when a» 
increases are discernable. The bottom left panel shows 
the evolution of the semimajor axis and how close the 
planet is to its Roche limit. As for WASP-4b, the thick 
dots indicate the end of the evolution where the peri- 
astron reaches the Roche limit. Despite the P^ 1 de- 
pendence of Q'„ (bottom right panel) and, therefore, the 
smoother plunging of the planet when the orbit has circu- 
larized (compared with the constant-Q* case), the eccen- 
tricity and the semimajor axis are decreasing extremely 
fast, especially after 1 Gyr. Thus, the age ranges when 
the radius of the planet, the eccentricity, and the semi- 
major axis simultaneously fit the measurements, are even 
narrower than for WASP-4b. Depicted by vertical orange 
segments, they are of the order of 50 Myr. Observing the 
planet in such a short interval of its life, questionable for 
WASP-4b, is even less likely for WASP-12b. Neverthe- 
less, as we described for WASP-4b in §4.2, using a 2a 
or 3cr criterion for fitting the observed properties of the 
system would similarly alleviate this problem. Finally, 



as in the case of WASP-4b, the radius when the orbit 
has circularized is clearly above the radius obtained with 
no tides. 

5. THE PLUNGING TIMESCALE UNCERTAINTY 

An intriguing issue about the transiting extrasolar gi- 
ant planets is their subsequent evolution. Tidal evolu- 
tionary equations show that their fate is tidal disruption 
as they inspiral into their host star (Rasio et al. 1996; 
Levrard et al. 2009; Jackson et al. 2009b, a; Ibgui & Bur- 
rows 2009; Miller et al. 2009). Once the orbit has cir- 
cularized, the plunging timescale r is described by an 
analytic expression, directly resulting from the evolution 
of a (eq. 6), where P is the orbital period: 

- (£)(£)'■ < 10 ' 

We define r as the remaining time for the planet to fall 
into its host star from the present time, assuming that 
its current eccentricity is zero. We thus choose in this 
section to, defined in §2.3, to be the present time, and 
therefore do to be the measured semimajor axis. 

This formula shows the sensitive dependence of r on 
the current semimajor axis of the system ao- The latter 
is, however, quite precisely determined by observations 
(see Table 1). At the same time, eq. (10) exhibits a 
linear dependence of r on <3'„, which is a very poorly 
constrained parameter. 

To illustrate our point, we plot on Fig. 3 the theoreti- 
cal evolution, starting from the present time, of the semi- 
major axes o(AU) of two transiting EGPs, HD 209458b 
and WASP-12b. The time (in Myr) is represented log- 
arithmically. Our aim is to emphasize the dependence 
of these evolution on Q'^ that follows the generic law de- 
scribed by eq. (7). The reference orbital period Pq is 
the observed period, Po = -Pmeasurcd- We explore three 
possibilities for 7, assumed to range from -1 to +1 (see 
§6): 1 (dotted curves), O(solid), -1 (dashed). We also ex- 
amine three possibilities for the (5 factor : 5 (red curves), 
6 (blue), 7(grecn). This enables us to roughly encom- 
pass the current observational and theoretical estimates 
of . The measured eccentricities of these systems are 
low enough to consider the tides raised on the planet to 
be negligible and, therefore, to ignore Q' p . The analytical 
formula for the evolution of a (eq. 6) is applicable and so 
is the above definition of the plunging timescale (eq. 10). 
We numerically checked this point by comparison with 
the integration of the full evolutionary equations of §2.2. 
The Roche limits of both systems are plotted; they mark 
the end of the evolutionary curves represented by thick 
dots. 

The linear dependence of t on Q'„ directly affects the 
evolutionary curves. For HD 209458b, the order of mag- 
nitude of the plunging timescale can be 0.5 Gyr (/3 = 5), 
5 Gyr (/3 = 6), or 50 Gyr (f3 = 7). The dependence on 
the 7 parameter is weaker, since the ratio -P/Pmeasured 
remains around 1. The trend is that a positive 7 re- 
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suits in the decrease of while the orbital period P 
of the planet decreases as it spirals into its host star. 
This results in an increase of the tidal torque exerted 
on the planet and, therefore, in accelerated plunging in 
comparison to the case with 7 = constant) and 
a fortiori to the case with 7=1 oc P^ 1 ). Includ- 
ing all these uncertainties, HD 209458b can plunge be- 
tween 0.5 and 60 Gyr from now, which is a 2-order-of- 
magnitude range. Note that Levrard et al. (2009) pro- 
vide the evolution curve of HD 209458b for one of the 
cases considered here: (3 = 6 and 7 = 1, (blue dotted 
curve). Their result is consistent with ours. By the same 
token, wc find that WASP-12b can plunge between 0.1 
and 100 Myr from now, which is a 3-order-of-magnitudc 
range. The much shorter timescale, in comparison with 
that for HD 209458b, is mainly due to the ratio of the 
semimajor axis (a = 0.0471 AU for HD 209458b and 
a = 0.0229 AU for WASP-12b, see Table 1) combined 

with the power dependence t oc aj 3 ^ 2 . Note that for 
WASP-12b, the relative influence of the 7 parameter 
compared with that for the (3 parameter is bigger than 
for HD 209458b. The evolutionary curves can overlap. 

This figure demonstrates that it is difficult to predict 
the evolution of the orbits of transiting EGPs, given the 
poor knowledge of the tidal dissipation factors in the host 
stars. However, the recently discovered WASP-18b (Hel- 
lier et al. 2009), which has an orbital period of only 
0.94 days, might experience a measurable tidal decay 
time. Hellier et al. (2009) state that if Q' t = 10 6 , the 
epoch of transit would shift by 28 s after 10 yr. There are 
several other effects that, in principle, could also cause 
changes in the timing properties of transits. Miralda- 
Escudc (2002) points out that torques due to both the 
quadrupole moment of the star and the gravitational per- 
turbations from a hypothetical companion planet could 
cause precession of the orbital plane and of the planet's 
periapse, with associated effects on transit time and du- 
ration. Furthermore, Rafikov (2009) suggests that both 
general relativistic apsidal precession and proper motion 
of the exoplanetary system with respect to our solar sys- 
tem could cause similar or greater changes in transit time 
and duration. 

6. CONCLUSIONS AND DISCUSSION 

We have presented in this paper some new general re- 
sults of the coupled radius-orbit evolutionary model de- 
scribed in Ibgui & Burrows (2009), and we have applied 
the model to the inflated planets WASP-4b and WASP- 
12b. We assumed a two-body gravitational and tidal in- 
teraction between the planet and its host star, coupling 
the planetary radius and the orbit evolution. We in- 
cluded the tides raised on the planet and the tides raised 
on the star. Stellar irradiation and a detailed planetary 
atmosphere are included. The fundamental result is the 
transient inflation of the planetry radius that temporarily 
interrupts its monotonic standard shrinking. An impor- 
tant point is that even though the current orbit of the 
planet has almost circularized, the radius of the planet 
can still be inflated due to an earlier episode of tidal heat- 
ing. This is why we stress that an inflated planet with 
an observed circular orbit can still have tidal heating as 
an explanation of its radius. Fixing the planet and star 



properties, the model is controlled by four free parame- 
ters, (Q' p , ej, ttj), that are the tidal dissipation factors 
in the planet and in the star, and the initial eccentricity 
and semimajor axis at the beginning of this two-body 
evolution. We stress the sensitive and nonlinear depen- 
dence of the evolutionary curves on these parameters. 

We have demonstrated that an increase of cither the 
core mass M coro , or Q' p , or a, results in a lower value of 
the radius inflation peak and in a delay of its appearance. 
The final semimajor axis is the same, whatever M COIC or 
Q' p , but is larger when a, is larger. At an earlier age, the 
planet with the larger core has the smaller radius, but 
this is opposite at later ages. 

We have enhanced our model by including an orbital 
period dependence of the tidal dissipation in the star, 
Q\ oc P 7 , — 1 < 7 ^ 1. Q' t drives the inspiral of the 
planet into its host star. 

Applications of our model to recently detected transit- 
ing inflated planets show that: 

• WASP-6b and WASP- 15b can be fit at solar opac- 
ity over Gyr age ranges. 

• We have not found an acceptable fit for TrES-4, 
at either solar, 3 x solar, or 10 x solar planet atmo- 
spheric opacity. 

• WASP-4b can be fit at solar opac- 
ity with, for example, the combination 
{Q' p ,ei,ai)=(10 8 - , 0.80, 0.050) and with 
Q'„ = 10 6 - 5 x (P/lOdays)- 1 . 

• WASP- 12b can be fit at solar opac- 
ity with, for example, the combination 
(Q^,e 4 ,a,)=(10 8 - , 0.73, 0.055) and with 

Q'„ = 10 6 - 5 x (P/lOdays)- 1 . 

For WASP-4b and WASP-12b, the ranges of ages that 
allow simultaneous fits of radius, semimajor axis, and 
orbital eccentricity, are very narrow, seeming to sug- 
gest that, if the two-body coupled evolutionary model 
described herein is in fact responsible for these plan- 
ets' inflated radii, then we are observing them at a spe- 
cial epoch in their evolution. However, relaxing the fit- 
criterion from lcr to 2cr or 3c would alleviate this appar- 
ent problem. 

Our results (in particular, for TrES-4) suggest that a 
coupled radius-orbit tidal evolution model might not on 
its own explain the radii of all the inflated transiting 
giant planets. An alternative scenario with stationary 
heating has been proposed (Ibgui et al. 2009) and ap- 
plied to all the planets discussed in this paper. Though 
not providing direct solutions to the inflated radii is- 
sue, this scenario constrains the ratio e 2 /Q' p for a given 
Mcorc- Finally, a combination of these two models could 
be imagined with a two-body interaction, followed by a 
quasi-steady low eccentricity phase due to perturbations 
by a second planet. 

The last point we make in this paper is the uncer- 
tainty of the plunging timescale during the spiraling of 
the planet into its host star. This timescale is strongly 
dependent on the semimajor axis; specifically, it depends 
on a to the 6.5 power. It also has a linear dependence on 
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Q'„ , which is a parameter that is uncertain by several or- 
ders of magnitude. We have shown that HD 209458b can 
plunge in between 0.5 and 60 Gyr from now, a 2-order-of- 
magnitude range, and that WASP-12b can plunge in be- 
tween 0.1 and 100 Myr from now, a 3-order-of-magnitude 
range. 

Ibgui & Burrows (2009) have suggested caveats to, and 
ways to improve, the model employed here. We close by 
noting several additional points. The orbital evolution 
equations depend on the theory of tidal dissipation inside 
gaseous planets and stars. Improvements to this theory 
might result in different evolutionary tracks (of R p , e, 
and a) from the ones presented in this paper. Further- 
more, we have noted that it is not strictly appropriate to 
apply these equations to model scenarios with large val- 
ues of orbital eccentricity. However, both for comparison 
with previous work and because the proper tidal theory 
remains unknown, our present approach is a valuable 
step in exploring the extent to which tidal dissipation 
might explain the radii of the inflated EGPs. Further ob- 
servations that might help to constrain this model and to 
discriminate between this and the stationary-state model 
of Ibgui et al. (2009) include both increasing the accu- 



racy of orbital eccentricity measurements and searching 
for companions to the transiting EGPs. These and other 
advances will help us progress toward a better under- 
standing of the puzzle of the inflated planets. 
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TABLE 1 
Planet and Orbital Data 



Planet 


a 

(AU) 


e 


P 

(days) 


M p 

(Mj) 


(M e ) 


Rp 
(Rj) 




Roche limit 
(AU) 


References 


HD 209458b 


0.04707t° :S™« 


< 0.028 


3.52475 


0.685l° ;°ll 


218l| 


1 qon+O- 024 
1 - JZU -0 .025 


1.000 


0.019 


1,2,3 


TrES-4 


U - UOiUO -0. 00167 


< 0.01 c 


3.55395 




2941^ 


1 7 oo+0.093 


2.371 


0.024 


4 


WASP-4b 


02340+ 00060 

U.UZO'±U_q 00060 


< 0.096 


1.33823 


1 007+0. 064 
L.ZOl _o.064 


a3 °-20 


1 Ofic+0.021 
l.ODO_ 021 


1.706 


0.015 


5,3 


WASP-6b 


421 + 0008 
u.u^zi_ 0013 


0.054iHJ 8 5 


3.36101 


S03+ 019 
u.ouo_ 03g 


160l? 2 


1 99A+ 051 
± -^^ ± -0.052 


0.463 


0.018 


6 


WASP- 12b 


0229+ 0008 


0.049tHl 5 5 


1.09142 


1 41+0.10 

1 -^ i -0.10 


448li 


1 7Q+0-09 
± - ' y -0.09 


9.098 


0.021 


7 


WASP- 15b 


04 99+ 0018 


0.0521HS 


3.75207 


542+ 050 

u - dq:Z -0.050 


17211* 


1 A98+ 077 
1 '* io -0 .077 


1.696 


0.022 


8 



References. — (1) Knutson et al. (2007), (2) Torres et al. (2008), (3) Madhusudhan & Winn (2009), (4) Sozzetti et al. 
(2009), (5) Winn et al. (2009a), (6) Gillon et al. (2009a), (7) Hebb et al. (2009), (8) West et al. (2009). 

a F p is the stellar flux at the planet's substellar point, in units of 10 9 ergs cm~ 2 s — 1 . 

6 According to Hartman (2009), the 68.3% confidence upper limit is 0.065 and the 95.4% confidence upper limit is 0.085. 
c For specificity, we adopt this value of 0.01. Knutson et al. (2009) derive a 3(7 upper limit |ecos(w)| of 0.0058, where lo is the argument 
of periastron. 



TABLE 2 
Host Star Data 



Star 


M* 
M 


R* 
R© 


(K) 


[Fe/Hl 
(dex) 


Age 
(Gyr) 


HD 209458 


1 ini+ 066 


1 19 r+0.020 
J - J " 4J -0.023 


6065122 


o.ool2;2 5 5 


q 1+O.8 


TrES-4 


1 404+ 066 

l.tV4_ 134 


1.846lH 9 8 6 7 


62001^ 


+o.i4i2;22 


2.91- 


WASP-4 


n Q9C+0-040 

U.»ZO_Q Q4Q 


0.912±Hi! 


5500±i°8 


-0 n3 +009 

u.uo_ 09 


6.511:2 


WASP-6 


0.88012:222 


0.87012:222 


54501 1 ™ 


-0.2012.22 


Ht 7 7 


WASP- 12 


1 QK + 0.14 

1 - od -0.14 


1 1-7+0.07 

1 - d '-0.07 


6300t»g 


+0.3012:?^ 


211 


WASP- 15 


1 1S+ ' 12 
1.18_ 12 


1 477+0.072 

' -0.072 


6300112° 


-o.i7l2:ll 


q q + 2.8 
°- 3 -1.3 



Notes — The references are the same as the ones used for the data in Table 1. The ages are less well constrained and should be taken 
with caution. 
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Fig. 1. — Examples of fitting (within 1<t of the measurements) evolutionary curves for WASP-4b at solar atmospheric opacity, assuming coupled evolution with tides of the planetary 
radius and orbit (ibgui & Burrows 2009). Top left, top right, and bottom left show (versus the age in Gyr) the co-evolution of the planetary radius R p (Rj), the orbital eccentricity 
e, and the semimajor axis a(AU). Moreover, bottom right are plotted logio(Q'„) (solid, units on left axis), where Q' t is the tidal dissipation factor for the star, and the orbital period 
P(days) (dashed, units on right axis). The pairs of orange vertical segments define the ranges of ages for which the fits are simultaneously obtained for R p , e, and a, whose measured 
values arc plotted along with the la limits. Notice how narrow these ranges are (< 0.2 Gyr) compared with the wide range of possibilities for the age of the system (~4.6 Gyr). Three 
cases are plotted. They differ in the initial semimajor axis a^= 0.050(blue), 0.052(green), 0.054(red). For the three of them, the initial eccentricity is ei=0.80, the tidal dissipation 
factor for the planet is a constant Q' p = 10 8 , and evolves as P -1 with Q'„ = 10 6 ' 5 at P = 10 days. Thick dots end the evolutionary curves at the age, for each case, when the 
periastron of the orbit p = a(l — e) reaches the Roche limit, drawn in bottom left. Bottom right, the dashed black horizontal curve indicates the measured orbital period of the system. 
See Section 4 for a discussion. 
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Fig. 2. — Same as Fig.l, but for WASP-12b. Note the very narrow range in age (~50 Myrs) for which the fits are simultaneously obtained for R p , e, and a. Note also how close the 
orbit of the planet is to the Roche limit, which is roughly 0.021 AU. Given the uncertainties in the determination of the semimajor axis a and the orbital eccentricity, the periastron 
p = a(l — e) of the orbit might be slightly smaller than the Roche limit. We find 0.0207 < p < 0.0229. See Section 4 for more discussion. 
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Fig. 3. — Semimajor axis a(AU) evolution of transiting EGPs driven by tides in the star. We assume tidal evolution of a at constant planet 
and star properties. The emphasis is on the dependence of this evolution on the tidal dissipation factor in the star Q' f , which is assumed 
to vary with the orbital period P, as Q' t = lO' 3 X (P/P m easured) 7 i where P me asured is the measured orbital period, (3 and 7 are a priori 
unknown factors quantifying the evolution of Q' t . We consider three possibilities for (3: 5(red), 6(blue), 7(green), and three possibilities for 
7: l(dotted), O(solid), -l(dashed). These choices enable us roughly to encompass current observational and theoretical estimates of Q' t . 
The effect of the tidal dissipation factor in the planet Q' p is negligible in these cases because of the low orbital eccentricities. We choose 
HD 209458b and WASP-12b, two planets evolving on disparate timescales. The planets end up falling into their parent stars and get 
disrupted at their respective Roche limits (indicated on the figure). However, the curves show the very large uncertainties in the timescale 
to reach this end. WASP-12b can plunge in between 0.1 and 100 Myr from now, a 3-order-of-magnitude range. HD 209458b can plunge in 
between 0.5 and 60 Gyr from now, a 2-order-of-magnitude range. The main source of uncertainty comes from the /3 factor. To summarize, 
this figure demonstrates how difficult it is to predict the evolution of the orbits of transiting EGPs, given the very poor knowledge of the 
tidal dissipation factors in the host stars. See Section 4 for more discussion. 



